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Abstract. We present XMM-Newton results on the temperature profiles of a volume-limited sample of galaxy 
clusters at redshifts z ~ 0.3, selected from the REFLEX survey (REFLEX-DXL sample). In the spectral analysis, 
where only the energies above 1 keV were considered, we obtained consistent results on the temperature derived 
from the pn, MOS1 and MOS2 data. Useful temperature measurements could be performed out to radii with 
overdensity 500 (rsoo) for all nine clusters. We discovered a diversity in the temperature gradients at the outer 
cluster radii with examples of both flat and strongly decreasing profiles. Using the total mass and the gas mass 
profiles for the cluster RXCJ0307. 0—2840 we demonstrate that the errors on the mass estimates for the REFLEX- 
DXL clusters are within 25% up to rsoo- 
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1. Introduction 

The number density of galaxy clusters probes the cosmic 
evolution of large-scale structure (LSS) and thus provides 
an effective test of cosmological models. It is sensitve to 
the matter density (f2 m ) and the amplitude of the cosmic 
power spectra on cluster scale (ag) (e.g. Schuecker et al. 
2003). Its evolution is sensitive to the dark energy (^a) 
(e.g. Vikhlinin et al. 2002). The most massive clusters are 
especially important in tracing LSS evolution since they 
are expected to show the largest evolutionary effects. In 
addition, the X-ray properties of the most massive clus- 
ters should be easier to describe in hierarchical model- 
ing since the structure of the X-ray emitting intracluster 
plasma is essentially determined by gravitational effects 
and shock heating. With decreasing cluster mass and in- 
tracluster medium (ICM) temperature, non-gravitational 
effects play an important role before and after the shock 
heating (Voit & Bryan 2001; Voit et al. 2002; Zhang & Wu 
2003; Ponman et al. 2003). Therefore, the most massive 
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clusters provide the cleanest results in comparing theory 
with observations. 



In this project we are analysing an almost 
volume-complete sample of thirteen X-ray luminous 
{Lx > 10 45 erg s _1 for 0.1 - 2.4 keV) clusters selected 
from the ROSAT-ESO Flux-Limited X-ray (REFLEX) 
galaxy cluster survey (Bohringer et al. 2001a) in the 
redshift interval z — 0.27 to 0.31. There is only a very 
small correction to the volume completeness with a 
well known selection function for Lx > 10 45 erg s _1 
at the higher redshift as described in Bohringer ct al. 
(2003; Paper I). With this REFLEX-DXL (Distant X-ray 
Luminous) sample we want to obtain a robust measure 
of the cluster abundance of this epoch, in particular to 
perform studies of the evolution of the cluster population 
by comparing these observations with more nearby and 
more distant clusters. A prime goal is to obtain reliable 
ICM temperatures of these clusters as a measure of the 
cluster masses (e.g. Evrard 1997). Since peculiarities in 
the cluster structure introduce a scatter in the mass- 
temperature relation and since in particular on-going 
cluster mergers can lead to a temporary increase in 
the cluster temperature and X-ray luminosity (Randall 
et al. 2002), we aim for a detailed study of the deep 
XMM-Newton observations described here. The clusters 
are also scheduled for a detailed spectroscopic study 
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Table 1. Compilation of some observational information on the nine REFLEX-DXL clusters. Col. (1): Cluster name. 
Cols. (2-3): Sky coordinates. Cols. (4-6): Net exposure time of MOS1, MOS2 and pn after cleaning for the flaring 
episodes. Cols. (7-9): Light curve cleaning upper limit. Col. (10): Hydrogen column density in units of 10 20 cm~ 2 
(Dickey & Lockman 1990). Col. (11): Revolution of XMM-Newton. Col. (12): Alternative name 



Cluster 




a(°) 


s n 


Exposure Time (s) 


Criteria (cts/lOOs) 


riu 


orbit 


Alternative 


(RXCJ) 




Eq. J2000.0 


MOS1 


MOS2 


pn 


MOS1 


MOS2 


pn 






name 


0014.3 - 


3022 


3.5837 


-30.3757 


15085 


15510 


10057 


23.4 


23.7 


61.2 


1.60 


270 


A2744 (AC118) 


0043.4 - 


2037 


10.8508 


-20.6225 


11253 


11248 


6318 


23.5 


25.0 


69.8 


1.54 


380 


A2813 


0232.2 - 


4420 


38.0717 


-44.3453 


11979 


11508 


7741 


22.2 


22.4 


63.4 


2.49 


474 




0307.0 - 


2840 


46.7667 


-28.6708 


12309 


12610 


8126 


21.9 


22.5 


56.8 


1.36 


218 


A3088 


0528.9 - 


3927 


82.2342 


-39.4636 


7097 


6806 


3297 


23.3 


23.2 


57.4 


2.13 


324 




0532.9 - 


3701 


83.2350 


-37.0260 


10374 


11191 


6527 


23.3 


25.3 


63.0 


2.90 


518 




0658.5 - 


5556 


104.5700 


-55.9600 


25339 


23365 


18307 


24.7 


23.6 


57.7 


6.53 


159 


1ES0657-558 


1131.9 - 


1955 


172.9858 


-19.9258 


11660 


11164 


8511 


22.3 


22.3 


57.8 


4.50 


286 


A1300 


2337.6 + 


0016 


354.4204 


0.2760 


12216 


11915 


7412 


22.4 


23.3 


69.1 


3.82 


365 


A2631 



of the cluster dynamics with the ESO-VLT-VIMOS 
instrument. 

The selection of the REFLEX-DXL sample and its 
properties are described in detail in Paper I. For all clus- 
ters in this sample the XMM-Newton observations have 
confirmed that the X-ray luminosity is dominated by 
diffuse thermal emission from the ICM of these sys- 
tems. Therefore, the REFLEX-DXL sample contributes 
a unique sample of X-ray luminous and consequently very 
massive clusters from roughly the same epoch, which are 
not only interesting as cosmological probes, but also for 
astrophysical studies like the statistics of cluster substruc- 
ture, galaxy evolution, Sunyaev-Zel'dovich observations 
and many other applications (e.g. Bohringer et al. 2001b). 

The estimate or derivation of the cluster mass is an 
essential step in almost all these studies. The mass can 
be either approximately estimated from the temperature 
(Evrard 1997), or determined from the temperature and 
density distributions of the ICM under the assumption 
of hydrostatic equilibrium of the intracluster gas (e.g. 
Cavalicrc & Fusco-Femiano 1976; Serio et al. 1981), or 
otherwise determined from the mass of the intracluster 
gas and the assumption of the universality of the cluster 
baryon fraction (e.g. White et al. 1993; Vikhlinin et al. 
2002). The first two methods require a robust determina- 
tion of the ICM temperature and a good understanding 
of the cluster structure for a reliable interpretation of the 
results. 

Therefore, it is the aim of this paper to establish 
a reliable method of spatially resolved temperature de- 
termination for the clusters in the REFLEX-DXL sam- 
ple and to derive temperature profiles for all the clus- 
ters. XMM-Newton with its superior sensitivity com- 
bined with its good spatial resolution provides the best 
means for such studies (Arnaud et al. 2002). Previously, 
large data sets on cluster temperature profiles have 
been compiled from ASCA (e.g. Markevitch et al. 1998; 
White 2000; Finoguenov et al. 2001a; Finoguenov et al. 



2002; Sanderson et al. 2003) and BeppoSAX observations 
(Molendi & De Grandi 1999; Ettori et al. 2002). 

This paper is structured as follows. In Sect. we de- 
scribe the background components, which are important 
to this study. Then we present a double background sub- 
traction method, which is developed to provide a precise 
background removal. In Sect.|3 we analyse the properties 
of the hot gas in the galaxy clusters, and show our analytic 
temperature model. Then we determine the total mass 
and gas mass fraction in the cluster RXCJ0307. 0—2840 
based on the precise temperature and gas density pro- 
files. In Sect. El we draw our conclusions. We adopt a flat 
ACDM cosmology with the density parameter fi m = 0.3 
and the Hubble constant Hq = 70 km s _1 Mpc -1 . Error 
bars correspond to the 68% confidence level, unless explic- 
itly stated otherwise. 

2. Method 

2.1. Data preparation 

Of the thirteen XMM-Newton observations of REFLEX- 
DXL clusters conducted so far, eight have sufficient qual- 
ity for the detailed studies described here. The remain- 
ing five clusters are heavily affected by soft proton flares. 
Some properties of these observations are described in 
Paper I. Re-observations of these targets have been al- 
located. An additional X-ray luminous REFLEX cluster 
RXCJ0307.0-2840 at z = 0.2578 was observed in this 
project and is also analysed here. It has very good ob- 
servational data and appears to be a very regular and 
symmetric cluster. We have therefore selected this object 
as an example to demonstrate our method of analysis. 

An overview on the observational data of the complete 
sample of thirteen plus one clusters is given in Paper I. 
In this paper we compiled further observational informa- 
tion on the sample targets in Tabled which includes the 
observational parameters of the data and the alternative 
names of these targets. 
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We use the XMMSAS v5.4 software for data reduction. 
The MOS and pn data were taken in standard Full Frame 
mode and Extended Full Frame mode, respectively. For 
all detectors, the thin filter has been used. 

Above 10 keV, there is little X-ray emission from the 
clusters due to the low telescope efficiency at these ener- 
gies; the particle background therefore completely domi- 
nates. The cluster emission is not variable, so any spec- 
tral range can be used for temporal variability studies of 
the background. Therefore, the 10-15 keV energy band 
(binned in 100 s intervals) was used to monitor the parti- 
cle background and to excise periods of high particle flux. 
In this screening process we use the settings FLAG = 
and PATTERN < 5 for pn, and PATTERN < 13 for 
MOS. As an example, Fig.Q]shows the 10-15 keV pn light 
curve of RXCJ0658. 5-5556. 

We reject those time intervals affected by flares in 
which the detector countrate (ctr) exceeds a threshold 
of 2cr above the average ctr, where the average and the 
variance have been interactively determined from the ctr 
histogram below the rejection threshold. A similar clean- 
ing criterion is applied to the screening of the background 
observation. We note, however, that the thresholds will 
be different for the source and background accumulations. 
Formal freezing of the cleaning criterion does not overcome 
the difference in the mean background ctr. In our analysis 
we searched through a number of background observations 
to find the one matching our target observations. The se- 
lection criterion is therefore to find the one with a similar 
cleaning threshold. 

As shown in Table Q the pn and MOS detectors all 
have their own similar cleaning thresholds for all obser- 
vations. In Sect. 12.21 we consider in detail how this back- 
ground behaviour affects the temperature determination. 

In the analysis of the pn data, we statistically remove 
the out-of-time effect by creating an out-of-time (OOT) 



Table 2. Parameters of the residual background models 
fitted in the 0.4-15 keV band. Col. (1): Cluster name. 
Cols. (2-4): Index of the "powerlaw/b" residual back- 
ground model for MOS1, MOS2 and pn. Cols. (5-7): 
Normalization at 1 keV of the "powerlaw/b" residual 
background model scaled to 1 arcmin 2 area for MOS1, 
MOS2 and pn in units of 10~ 4 cts s" 1 keV^ 1 arcmin . 



Cluster 




Index 




Normalization 


(RXCJ) 


MOS1 MOS2 pn 


MOS1 MOS2 pn 


0014.3-3022 


1.47 


1.48 


1.95 


1.21 


1.73 


6.49 


0043.4-2037 


1.71 


1.43 


1.61 


1.57 


1.28 


8.50 


0232.2-4420 


1.26 


1.52 


1.54 


1.26 


1.73 


6.60 


0307.0-2840 


1.28 


1.43 


1.82 


0.89 


1.21 


2.80 


0528.9-3927 


0.80 


0.96 


1.56 


0.93 


1.12 


3.36 


0532.9-3701 


1.60 


1.67 


2.08 


0.32 


1.51 


3.89 


0658.5-5556 


1.52 


1.64 


1.95 


7.66 


3.99 


15.26 


1131.9-1955 


1.98 


2.34 


3.19 


1.73 


1.89 


5.96 


2337.6+0016 


1.24 


1.47 


1.43 


1.71 


1.58 


8.74 



event list file and using the XMMSAS products such as 
images and spectra to subtract it. The present observa- 
tions have been taken in the Extended Full Frame mode 
(frame time=199 ms). At this mode, the fraction of 
OOT effect amounts 2.32%, which we used to nor- 
malize the XMMSAS OOT product before we sub- 
tract it from the XMMSAS normal product. 

2.2. Background analysis 

The purpose of the background analysis described in the 
following is to find a suitable "blank field" observation 
to be used for the background subtraction and to further 
analyse the difference between the target and background 
to take this residual background into account. 

The background has several components which exhibit 
different spectral and temporal characteristics. In the low 
energy band (< 0.3 keV), the instrumental background is 
dominated by electronic noise consisting of a large number 
of small amplitude events added up during each frame ac- 
cumulation (Read 2002). This noise depends on the read- 
out frequency of the cameras and is sensitive to the en- 
ergy offset of each individual pixel. Energetic particles 
produce several line and continuum components in the 
background, which can be further subdivided into time 
variable and constant components. The constant compo- 
nent has been intensively studied by Lumb et al. (2002) 
and Read & Ponman (2003). This component can be re- 
moved using the so called blank field observations. De 
Luca & Molendi (2001) report some evidence for a sec- 
ular evolution of the background level on a half year time 
scale. Therefore, background accumulations close to the 
date of the target observation are more suitable. In addi- 
tion, variations in the instrumental background on a much 
shorter time scale have been seen (e.g. De Luca & Molendi 
2001). Part of such periods with an increased background 




8.855X10 7 8-856X10 7 8.857X10 7 8.858X10 7 

time (s) 



Fig. 1. pn light curve of RXCJ0658. 5-5556 in the 10-15 keV 
energy band. Time is measured in second relative to the 
XMM-Newton internal clock. 
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arc rejected through the analysis of the light curve (e.g. 
Read & Ponman 2003). However, we sometimes still ob- 
serve a residual enhancement of the background associ- 
ated with an increase in the quiet flux of soft protons. The 
typical time-scale for the variation of this 'quiet' compo- 
nent is comparable to or exceeds the typical observational 
time-scale. Therefore, it is in general not possible to re- 
move observational intervals affected by this enhancement. 
Chen et al. (2003) describe an example of such an obser- 
vation where quite different 'quiet' background levels are 
observed before and after a flare, respectively. 

The photon background consists of foreground emis- 
sion from the Galaxy as well as the Cosmic X-ray 
Background (CXB). Observations of the blank field also 
contain both components, provided the accumulations are 
done with the same instrumental set-up (e.g. with a par- 
ticular filter) and the spectra of the X-ray background are 
the same for both the target and the blank field. This is 
only guaranteed for the CXB and the emission from the 
Galaxy halo. Since the Galaxy also plays a role as absorber 
and foreground emitter, it is important to have a similar 
absorbing column density for both the target and back- 
ground accumulations. Additionally, there are some extra 
Galactic components, that display spatial variations. To 
constrain them, one has to look into the ROSAT All-Sky 
Survey data around both the background accumulations 
and the target. Observations with normal conditions of 
the Galactic emission are referred to the quiet Galactic 
zones. So are most of the X-ray background accumula- 
tions (Lumb et al. 2002; Read & Ponman 2003). Some 
source removal is performed on the existing background 
accumulations. This changes the shape and intensity of 
the residual CXB. Therefore, a similar source removal has 
to be performed in the analysis of the target. 

Several available XMM-Newton pointings have been 
investigated. We conclude that the XMM-Newton point- 
ings in the Chandra Deep Field South (CDFS) have sim- 
ilar background conditions as most targets and sufficient 
exposure time. Therefore, the CDFS is a good candidate 
for the background for our study. However, there is still 
a small difference between the background of our sources 
and the CDFS, e.g. the background ctr in the target ob- 
servations is slightly higher (10—20%) than in the CDFS, 
which we ascribe mostly to the contamination by soft pro- 
tons. 

We have carefully planned this cluster study such that 
the radii in which the cluster emission can be observed 
extent up to spherical overdensities of 500, i.e., the ratio 
of the mean density of the dark halo with respect to the 
redshift-dependent critical density p c tn,{z). This is the re- 
gion to which the cluster X-ray emission is expected to be 
essentially confined, covering about half of the field of view 
(FOV) of the XMM-Newton telescope. This enables us to 
extract a source spectrum from the background region of 
the target field for comparison with a background spec- 
trum from the background field extracted with the same 
detector coordinates. Our residual background modeling 
procedure consists of analysing such regions not affected 



by cluster emission. We assume little or no vignetting of 
the soft proton induced background, as suggested by re- 
cent studies (Lumb et al. 2002). Spectra are extracted 
from the 9.2-11.5' region from the pointing center for our 
source observations and background candidate observa- 
tions. In a first step, we compare the spectra extracted 
from the outer regions in the 0.4-15 keV band between 
the sources and the background pointings. The residual 
background signal found after the subtraction of the back- 
ground field from the target field in this outer area is then 
modeled by a power law spectrum (model "powerlaw/b" 
in XSPEC, a power law background model which is con- 
volved with the instrumental redistribution matrix but not 
with the effective area). We use this model to account 
for the excess soft proton background in our observations 
as compared to the background field. This residual power 
law background model is introduced over the whole energy 
range ("wabs*mekal+powerlaw/b" , an emission spectrum 
from hot diffuse gas, e.g. Mewe et al. 1985, considering 
the Galactic absorption and modeling the residual back- 
ground by a power law), which yields the correct shape 
of the background component after the combination with 
the blank field background (double background subtrac- 
tion method). During this procedure, the normalization 
of the residual background component is always scaled to 
the area of the source extraction region. 

In Fig. |3 we present examples of the residual back- 
ground. The parameters of the residual background mod- 
els fitted in the 0.4-15 keV band for the clusters are given 
in Table El The residual background in pn is higher than 
in MOS because the larger thickness of the pn-pixels leads 
to a higher sensitivity to the particle flux. Since the sub- 
traction of the residual background is only a second or- 
der correction to the data and because of the large noise 
in the residual background data, we are not attempting 
a perfect model fitting, but approximating the data by 
a simple power law model. The uncertainty in the nor- 
malization is, anyway, within 5% and 10% for MOS and 
pn, respectively. The correction due to the residual back- 
ground makes only a 1-4% effect in the cluster signals and 
an up to 10% effect to the temperature determinations for 
cluster radii r < 4'. At larger radii the effect of the correct 
background effect is large as shown in Fig. [3] but the un- 
certainty in the approximation of the residual background 
- a third order effect - is still small. 

To recover the correct spectral shape and normaliza- 
tion of the cluster component, we need both the response 
matrix file (rmf) and auxiliary response file (arf). The 
following need to be taken into account in either 
rmf or arf: (i) Pure redistribution matrix giving 
the probability that a photon of energy E, once 
detected, will be measured in data channel PI. (ii) 
Quantum efficiency (with closed filter position) of 
the CCD detector, (iii) Filter transmission, (iv) 
Geometric factors such as bad pixel corrections 
and gap corrections (around 4%). (v) Telescope 
effective area as a function of photon energy, (vi) 
Vignetting correction to effective area for off-axis 
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pointings. We choose the rmf which corresponds 
to (i) and (ii) (with closed filter position). The pn 
rmf that corresponds to this choice has a '.closed' key- 
word in their naming conventions. For the MOS detec- 
tors we use software to generate such files, kindly pro- 
vided by S. Sembay. The arf corresponds to (iii), 
(iv) and (v) and is made according to the aver- 
age flux detected in the different extraction an- 
nuli, which takes (vi) into account. It is created us- 
ing the XMMSAS based program 'clarf ' by A. Finoguenov. 
Furthermore, in our analysis we apply the logarith- 
mically spaced radial bins, which provide nearly 
the same flux per bin. The importance of the 
scaterring due to the point-spread-function (PSF) 
is therefore suppressed. Using the XMM-Newton 
PSF calibrations by Ghizzardi (2001) we have es- 
timated the loss fraction of the flux: 20% for the 
central bins (0.5') and less than 10% for the other 
bins (> 1') with energy dependent effects being 
negligible. 

In summary the spectral analysis is performed in two 
steps, (i) A model for the residual background (back- 
ground difference) is obtained in XSPEC from a compari- 
son of the outer region of the target and background fields 
(see Fig. Eland Table |2J). (ii) The spectral modeling is per- 
formed in XSPEC with the cluster region as source data, 
the CDFS as background and the residual background as 
a second model component with model parameters fixed 
to the values found in step (i) (see Fig. [5j. 

In Fig. |3| we present the temperature profiles resulting 
from this background subtraction method fitted in the 0.4- 
10 kcV band for the example object, RXCJ0307.0-2840. 
We actually compare the results of the two step back- 
ground subtraction considering the residual background 
with the simple one step background subtraction. In Fig. 01 
we provide the comparisons of the results from MOS1, 
MOS2 and pn fitted in the 0.4-10 keV band. One notes 
that the residual background subtraction provides results 
in which all three instruments tend to show a slightly bet- 
ter consistency. The upper limit of the uncertainties be- 
tween the instruments goes down from 90% to 15% after 
the residual background subtraction. A detailed treatment 
of the background does not completely remove the differ- 
ences between the instruments. We found systematically 
lower temperatures obtained with pn compared to MOS1 
and MOS2 (partially). Because the pn detector is sensi- 
tive to the soft component, the pn measurements are easily 
affected by the soft band, which results in the lower tem- 
peratures given by pn compared to MOS. Since the tem- 
perature estimates of the pn are more strongly dependent 
on the soft energy band compared to MOS, we have car- 
ried out a spectral analysis in the 1-10 keV band. Despite 
the larger error bars, all temperature determinations in 
the central three bins of the cluster become higher, once 
the 0.4-1 keV band is excluded. In the following we will 
systematically investigate the effect of an energy band se- 
lection in the temperature and mass estimate. 



3. Results 

3.1. Redshift, mean temperature, and metallicity 

In a first step of the data analysis we derive global prop- 
erties of all nine galaxy clusters with good XMM-Newton 
data. Since a fraction of the clusters have cooling cores, 
dense central regions with lower temperature and so-called 
cooling flows, while the others do not display this phe- 
nomenon, we derive global temperatures including and 
excluding these regions. In addition the signal-to-noise de- 
creases fast in the outer regions of the clusters. Therefore, 
the global temperature was determined in the r < 8' re- 
gion and alternatively in the 0.5 < r < 4' region (see 
Fig. [5J . The global temperatures determined in both re- 
gions show some differences. The explanation is partially 
revealed by the temperature profiles. The metallicities in 
both zones are very similar. 

After subtracting the background and applying the 
rmf and arf, we fit the spectra in XSPEC using the 
one and two step correction models ( "mekal*wabs" , 
"mekal*wabs-(-powerlaw/b"). The fit using the latter 
model is better. For the regions covering radii of 0.5 < r < 
4' and r < 8', respectively, we use the 0.4-10 keV and 1- 
10 keV energy bands. Furthermore, we exclude the regions 
of substructure and/or some very bright point sources for 
several clusters throughout the procedure (cf . Table OJ) . 
The temperature of the small and large substructures in 
the well-known post merger cluster RXCJ0658.5— 5556 are 
8.3j^4 keV and 15.0l^g keV from the double background 
subtraction method in the 2-12 keV band. We exclude 
the small substructure to measure the global temperature. 
The flux within the region we excluded contributes 60% 
(10") to 85% (30") to the total flux of the point source. 
The temperature measurements vary within 10% after 
the subtraction. 

The redshifts obtained from the X-ray data (see 
Table which are ascribed to the ICM, are consistent 
with the redshifts obtained from optical spectra of indi- 
vidual cluster galaxies (see Bohringer et al. 2003) except 
for RXCJ2337. 6+0016. The optical redshifts contain a he- 
liocentric correction, while the uncertainty in the X-ray 
determined redshifts is one to two magnitude larger than 
this correction and thus no correction was made. We be- 
lieve that the optical redshift of RXCJ2337.6+0016 with 
5 coincident cluster galaxy redshifts is more reliable and 
accurate than the X-ray result at this stage. We plan to 
obtain further information on this object to resolve this 
discrepancy. 

The measurements of the global temperatures are sum- 
marized in Table 0] Similar to the measurements for the 
0.5 < r < 4' region, the measurements for the r < 8' 
region in the 0.4-10 keV band are lower than in the 1- 
10 keV band. Therefore, we only provide the comparison 
of the measurements for the 0.5 < r < 4' region fitted in 
two bands. The results obtained for the full and restricted 
spatial zones are consistent within 1-2ct (formal errors) 
within 15% uncertainties. 
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The global temperatures obtained from model fits to 
the larger spectral range 0.4-10 keV are always lower 
compared to the temperatures obtained from 1-10 keV. 
To check if the discrepancy is partially due to residual 
Galactic emission, we undertake the following test. We 
extract the spectra from the inner (hereafter "A") and 
outer (hereafter "B") parts of the background region in 
both the background (hereafter "bkg") and target (here- 
after "src") observations. If there is some difference in 
the Galactic emission between the background and target 
observations, there must be some residual Galactic emis- 
sion after we subtract the spectrum "B(src) — B(bkg)", 
scaled by the area of the region A, from the spectrum 
ll A(src) — A(bkg)" because of the vignetting effect of the 
Galactic emission. We apply a thermal emission spectral 
model ("apec ") with all parameters free in XSPEC to 
fit this residual emission. We found that the temperature 
of this component is around 0.2 keV and the redshift is 
close to zero. In the following analysis, we fix the tempera- 
ture to 0.2 keV, the abundance to solar abundance and the 
redshift to zero, and obtain the normalization of this com- 
ponent. We re-analyse the global properties of the clusters 
introducing the residual Galactic emission normalized by 
the area and vignetting effect. Since the difference in the 
temperatures determined in the 0.4-10 keV and 1-10 keV 
bands still remains, we report the temperature measure- 
ments setting the normalization of the residual Galactic 
emission to its upper limit (cf. Table 0J|. However, this 
component only makes a < 1 % effect in the cluster sig- 
nals. It is thus clear that the discrepancy is not or not 
mainly due to the residual Galactic emission. This dis- 
crepancy will again be discussed below. 



Table 3. Regions with substructures and point sources 
excluded from the analysis. Col. (1): Cluster name. Col. 
(2-3): Center of the circle in sky coordinates for J2000.0. 
Col. (4): Radii. Col. (5): "Yes" if there is an optical point- 
like counterpart in a digitized optical survey (e.g. DSS2). 



Cluster (RXCJ) 


a(°) 


S (") 


r (") 


opt 


0014.3-3022 


3.6306 


-30.3754 


15 


Yes 




3.5198 


-30.4154 


15 


Yes 




3.5172 


-30.4176 


10 


Yes 




3.5353 


-30.3654 


30 




0232.2-4420 


38.1566 


-44.3634 


20 


Yes 


0528.9-3927 


82.2427 


-39.4647 


15 


Yes 




82.3071 


-39.4809 


15 


Yes 




82.1450 


-39.3867 


15 


Yes 




82.1435 


-39.3724 


10 


Yes 




82.1134 


-39.3713 


10 


Yes 


0532.9-3701 


83.1957 


-36.9415 


15 






83.1506 


-37.0285 


10 




0658.5-5556 


104.5884 


-55.9413 


20 





3.1.1. Comparison with previous results 

Mushotzky & Scharf (1997) have measured for 
RXCJ0014.3-3022 a temperature of 12.Q8tJ;fjj keV 
(2cr errors) with ASCA data. Horner (2001) 
presents the temperature for the same cluster of 
9.61±g;f| keV, for RXCJ0232. 2-4420 of 7.19±EJ;3 2 ,keV, 
for RXCJ0307.0-2840 of 6.71 j^keV, and for 
RXCJ1131.9-1955 of 8.26±g;|l keV, all with 2a er- 
rors based on ASCA data. Lemonon et al. (1997) 
measured a temperature for the post-merger cluster 
RXGJ1131.9-1955 of 5 ± 3 keV using ROSAT PSPC 
data. They also found some evidence for a temperature 
gradient. 

For the cluster RXCJ0658.5— 5556 the published tem- 
perature measurements are less consistent. Tucker et al. 

(1998) measured a temperature of 17.4 ± 2.5 keV (2a 
error ) with ASCA. Andreani et al. (1999) obtained 
14.51?;? keV using both ASCA and ROSAT data. Yaqoob 

(1999) measured 11 ~ 12 keV with ASCA data. He found 
that a temperature of ~ 17 keV can be artificially obtained 
if the true spectrum has a stronger low-energy cut-off than 
that for Galactic absorption only. 

We notice that the ASCA spectra of 
RXCJ0658.5-5556 have in the 0.5-1 keV band only 
a few data points with large error bars. The differences 
in the temperature measurements described above come 
therefore clearly from the inclusion or exclusion of this 
soft part of the spectrum. Studies of nearby clusters 
suggest that putative non-thermal and warm thermal 
components are important at softer energies, while for 
rich clusters like Coma, the ICM dominates the X-ray 
emission up to 25 keV (Fusco-Femiano et al. 1999). 
Non-thermal emission dominates the emission at energies 
above 3 keV only in some of the groups of galaxies 
(Fukazawa et al. 2000). Since the expected temperatures 
of rich clusters are higher than 4 keV, we consider the 
temperature determination in the hard energy band 
as a more robust measure of the dominant gas mass 
component of a cluster, which traces the total mass. For 
RXCJ0658.5— 5556 we decided to restrict our tempera- 
ture determinations to the energy range 2-12 keV (see 
Tables H and 0. 

Our results on RXCJ0658.5— 5556 are consistent with 
Chandra measurements obtained by Markevitch et al. 
(2002) yielding a temperature of 14.8j^'2 keV from a fit 
of a spectrum extracted from the central r < 3' region. 
Although the authors show that U-qT > 15 keV in some 
parts of this cluster, the temperatures have quite large er- 
ror bars of 7 keV so that we regard this finding as not very 
significant. 

Additionally, Markevitch et al. (2002) fixed the value 
of the galactic hydrogen column density to nu = 4.6 x 
10 20 cm -2 which is significantly lower than jih = 6.5 x 
10 20 cm~ 2 obtained from Dickey & Lockman (1999). To 
check this result with our XMM-Newton data we set all 
parameters to be free to fit the spectrum extracted from 
the annulus region covering radii of 0.5 < r < 4' in the 
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0.4-12 keV band. In this case a high temperature of k^,T = 
18.8 ± 2.1 keV was obtained while the jih went down to 
an unrealistically small value of 1.8 ± 0.5 x 10 20 cm -2 , al- 
though the remaining parameters were still relatively rea- 
sonable. Therefore, we decide to fix nn = 6.5 x 10 20 cm~ 2 , 
but exclude the soft band (0.4-2 keV). In this case we ob- 
tained stable results for the temperature, metallicity, and 
redshift (see Table UJ. No significant metallicity gradients 
were found in our analysis. 

We notice that the temperature of RXCJ0528. 9-3927 
(also in other clusters) changes significantly with the low 
cut-off of the energy band used in the fit. In Fig. 
we thus used the X-ray spectra to test the energy band 
dependence and possible method dependencies by com- 
paring the temperature measurements versus low energy 
band (low-E) cut-off from two different methods (the dou- 
ble background subtraction method and the method ap- 
plied in Arnaud et al. 2002, i.e. we use the standard 
XMMSAS command 'evigweight' to correct vignetting, 
use 'arfgen' and 'rmfgen' to create on-axis arf and rmf, 
and apply the blank sky provided by Lumb 2002). For 
comparison, we applied the following models to fit the 
spectra in the 0.4-10 keV band after double background 
subtraction. We obtain /cbT = 7.69lg 46 keV and njj ~ 
— 0.1 x 10 20 cm~ 2 using a single-phase temperature 
model ("wabs*mekal-|-powerlaw/b") with free nn. We 
obtain k-sT = 10.12 ± 3.13 keV using two component 
thermal model ("wabs*(mekal+powerlaw)+powerlaw/b") 
with fixed n H - We also obtain fc B ?i = 9.34+^q4 keV and 
&bT2 ~ 0.49 — 1.34 keV using two thermal component 
model ("wabs*(mekal+apec)+powerlaw/b"), in which we 
fix the redshift of the soft component to the redshift of the 
cluster. The metallicity and redshift measurements among 
the different modelings and different low-E cut-off vary 
within 5%. The results presented in Fig. suggest some 
influence of the low energy band on the temperature mea- 
surements. As we discussed above, the results obtained in 
the harder energy band should recover the correct cluster 
temperature. Similar phenomena are also found for A1413 
(Pratt & Arnaud 2002) using XMM-Newton data, and 
Coma, A1795 and A3112 (Nevalainen et al. 2003) based 
on the comparison of XMM-Newton and ROSAT PSPC 
observations. Nevalainen et al. interprete this as a soft ex- 
cess, possibly due to a 'warm-hot' intergalactic medium. 
We will analyse this feature of our sample in more detail 
in a forthcoming paper. 

3.2. Temperature profiles 

We have already noted that the differences between the 
global temperatures of the regions covering radii of 0.5 < 
r < 4' and r < 8', respectively, in Table 0] are possi- 
bly caused by systematic temperature gradients. For a 
more detailed study of the temperature profiles we de- 
vide the cluster regions into the five annuli 0-0.5', 0.5-1', 
1-2', 2-4', and 4-8' (cf. Fig. - Note that in the spectra 
extracted from the outermost rings of RXCJ0014.3— 3022 



and RXCJ0528.9— 3927 we ignore a narrow energy band 
containing the residual background around the 1.49 keV 
Al line (Freyberg et al. 2002). 

Table shows the temperature profile catalogue 
of the clusters from the spectral analysis fitted by 
one or two step background subtraction ( "mekal*wabs" , 
"mekal*wabs+powerlaw/b"). We use both the > 0.4 keV 
and > 1 keV energy bands except for RXCJ0658.5— 5556. 
We apply the 2-12 keV band for this high temperature 
cluster. The temperature profiles of the clusters are also 
shown in Fig. |SJ 

For most of the objects, the temperature can be mea- 
sured out to rsoo ■ The overall temperature profile is char- 
acterized by a rather moderate decrease towards the cen- 
ter, and a decrease towards the outer regions, yet on dif- 
ferent levels, including no decrease at all for some of the 
clusters. This confirms a suggestion of Finoguenov et al. 
(2001b) that the differences in the behaviour of the tem- 
perature profiles in the outskirts of clusters have a statis- 
tical origin, rather than simple reflections of measurement 
errors. 

To test the validity of the results without a geometri- 
cal deprojection in our analysis, we apply the deprojection 
model provided in XSPEC ("projet") to study the de- 
projection effect for RXCJ0307.0-2840. This model per- 
forms a three dimension to two dimension projection of 
shells onto annuli. It is assumed that the inner boundary 
is specified by the outer boundary of the next annulus 
in. In the "projet" model, for each shell in a combined 
fit to all annuli spectra simutaneously, the contribution of 
each ellipsoidal shell to each annulus is determined and the 
spectral fitting results are then determined. In this fitting 
the outer shells are not affected by the emission from the 
inner shells. Similar work has been described by Pizzolato 
et al. (2003). Figure |H1 presents the temperature profiles 
from the spectral fit with and without projected modeling 
in the 0.4-10 keV energy band. The temperature gradient 
becomes slightly more significant when the geometrical de- 
projection effect is taken into account. The differences are, 
however, within the error bars and can thus be neglected. 
The relatively small effect of the deprojection is due to the 
steep surface brightness profiles of clusters which strongly 
reduce the influence of the emission from the outer shells 
on the observed spectra of the central regions. Therefore, 
the application of deprojection gives really significant im- 
provement only if the count statistics is very high (e.g. 
Matsushita et al. 2002). 

Systematic differences in the temperature profiles 
caused by the inclusion of the 0.4-1 keV energy band 
in the spectral analysis are not the same among 
the clusters. Cluster RXCJ1131.9-1955 is not af- 
fected at all, RXCJ0014.3-3022, RXCJ0307.0-2840 and 
RXCJ0528. 9-3927 are affected in the center, while 
RXCJ0043.4-2037 and RXCJ0232.2-4420 are affected in 
the outskirts. Since the instrumental setup used to ob- 
serve this sample is the same, it hardly is an instrumental 
artifact. However, more detailed analyses are needed in 
order to distinguish between the Galactic and extragalac- 
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Fig. 6. From left to right and from top to bottom, temperature profiles of RXCJ0014.3— 3022, 
RXCJ0043.4-2037, RXCJ0232.2-4420, RXCJ0307.0-2840, RXCJ0528.9-3927, RXCJ0532.9-3701, 
RXCJ0658.5-5556, RXCJ1131.9-1955 and RXCJ2337.6+0016. All the clusters except RXCJ0658.5-5556 
are fitted in the 0.4-10 keV band with (dashed lines) and without (dotted lines) the residual background 
subtraction. And the solid lines present the case with the residual background subtraction but fitted 
in the 1-10 keV band. The corresponding lines show the temperature profiles. Temperature profiles of 
RXCJ0658.5— 5556 are fitted in 2-12 keV with (solid lines) and without (dotted lines) the residual back- 
ground subtraction. 



tic origin of this component in the outskirts of some of the 
clusters in our sample (e.g. Finoguenov et al. 2003). 



3.3. Modeling RXCJ0307. 0-2840 

We use RXCJ0307.0-2840 as an illustrative example to 
demonstrate the accuracy of measurements of the total 
gravitating cluster mass and the gas mass fraction attain- 
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Table 4. Global temperatures, metallicities and redshifts of REFLEX-DXL clusters. Col. (1): Cluster name. Col. 
(2): Radius of annulus in arcmin. Col. (3): Energy band for fit. Col. (4): X-ray temperature measurements. Col. (5): 
Metallicity in solar abundance. Col. (6): Redshift obtained from the X-ray spectrum. Col. (7): \ 2 P er degree of freedom 
(d.o.f.). Col. (8): Redshift obtained from optical spectra as given in the REFLEX catalogue (see Bohringer et al. 2003). 
All X-ray spectra are fitted by the XSPEC model "mekal*wabs±powerlaw/b" . 



Cluster 


Region 


iinergy band (kevj 


kb.1 (kev) 


Z (Z ) 


— ray 


X /d.o.f. 


■^opt 


KAO J 0014. o-o022 


< r < 8 


i in 
1 — 1U 


q re+U.43 
O-OO-0.29 


A O /I 1 A A r 

0.24 ± 0.05 


A OA A _l_ A AAO 

0.294 ± O.OOo 


OAO A /OC\T 

oy2.9/o9/ 


A OACC 

O.oOdd 




0.5 < r < 4' 


0.4 — 10 


7 r-l +0.20 
' - oi -0.20 


0.22 zb 0.04 


A OTC 1 A AAA 

0.276 ± 0.009 


OOO O /0^71 

383.2/371 




* 






7 RQ+0-21 
' ■ DO -0.20 


A OO 1 A A A 

0.22 ± 0.04 


A 17/1 1 A AAT 

0.274 ± 0.007 


OTT O /OTA 

377.2/379 








1 — 10 


on+0.43 
8 - /y -0.32 


A OO 1 A AC 

0.22 ± 0.05 


A A T/" 1 1 A A 1 "1 

0.276 ± 0.011 


a a a /oro 

269.4/258 




KAGJ004o.4-20o7 


< r < 8 


1 — 10 


7 CA+ ' 47 

1 -OU_ 40 


A OO 1 A AC 

0.23 ± O.Oo 


A OAA 1 A AT ET 

0.309 ± 0.015 


1AO A /OOO 

198.0/238 


A OAO A 

0.2924 




0.5 < r < 4' 


0.4 — 10 


K qq+0.22 
D.88_ 2i 


A O A 1 A A/ 1 

0.24 ± 0.06 


A OAO 1 A A A/" 1 

0.303 ± 0.006 


1 AP7 O /1AT 

197.8/197 




* 






E ac+0.23 


0.25 ± 0.06 


A OAA 1 A A 1 1 

0.302 ± 0.011 


1A1 1 /lA^T 

191.1/197 








1 — 10 


ft CI +0.43 


A OO _l_ A AT 

0.23 ± 0.07 


A OAA 1 A Al O 

0.300 ± 0.012 


1 O A A 1 1 OA 

124.4/129 




nvri Tnoon o a ion 

HXC J 0232.2-4420 


< r < 8 


1 — 10 


ft 7A+0.24 

6-70I 025 


A OO 1 A A f 

0.32 ± 0.05 


A ri7T 1 A AAA 

0.275 ± 0.002 


AAA A /OOO 

292.0/332 


A A O O /"' 

0.2836 




a e ^ „ ^ ai 
0.0 < r < 4 


A A 1 A 
U.4 — 1U 


a qq+0.22 
°- oo -0.17 


A OO _J_ A AE 

U.22 ± U.05 


A OAE _l_ A AAO 

0.295 ± O.OOo 


OO/I 0/007 

384.2/28/ 




* 






a O7+0.20 

6.27I 017 


A OO 1 A A f 

0.23 ± 0.05 


A A A/" 1 1 A A"l A 

0.296 ± 0.010 


o i— 'O f /OOT 

372.6/287 








1 — 10 


7 co+0-40 
'■ D/ -0.33 


A O A 1 A A f 

0.24 ± 0.05 


A A A/" 1 1 A AAA 

0.296 ± 0.009 


1AP* A /lAr" 

195.4/195 




n vr 1 TAOH7 a no/in 
KAOJ0o0j\0-2o40 


a ^ ~. ^ q/ 
< r < 8 


1 1 A 

1 — 10 


a -I 7+0.24 

6.17_ 27 


A OO _1_ A AC 

0.28 ± O.Oo 


A O/ll _L A AAO 

0.241 zh 0.002 


OTA O /OOI 

2 /O.o/2ol 


A OE 70 

0.25/8 




0.5 < r < 4' 


0.4 — 10 


ft 1n +0.20 

6.10_ 022 


A OO 1 A AC 

0.28 ± 0.06 


A O A A 1 A AAO 

0.244 zb O.OOo 


275.0/251 




* 






ft i 7+0.21 

6.17I 021 


A OO 1 A AC 

0.28 ± 0.06 


A O A A 1 A AA A 

0.244 zb 0.004 


r\ry rr 1 /o 1 

255.1/251 








1 — 10 


ft ftq+0.34 

6.63I 031 


A OA 1 A A^7 

0.29 ± 0.07 


A O A O 1 A AA A 

0.24o zb 0.004 


1 7C C 1 1 C7 

176.5/167 




DVPTHCOO A OA07 

kaojod2o. y-^y2 r 


a ^ ~. ^ q/ 

< r < a 


1 1 A 

1 — 10 


C A7+0.59 


A 1 _J_ A AO 

O.ol ± O.Oo 


A OCO _l_ A Al O 

0.2t)2 zb 0.012 


IftO O/IEA 

163.2/109 


A OOOA 

0.2839 




a E ^ „ ^ a/ 
0.0 < r < 4 


A A 1 A 

U.4 — 1U 


ft ^n+0.29 

6.49_ 27 


A OA _1_ A AO 

O.oO ± O.Oo 


A O0 1 _l_ A Al 1 

0.2ol zb 0.011 


IftOO/l/IO 

162.8/148 




* 






« ,(0+0.28 
6 - 42 -0.27 


A 1 J A AO 

0.31 ± 0.08 


A OOI 1 A ATA 

0.2ol zb 0.010 


A r. A O / A A O 

154.8/ 148 








1 — 10 


t cc+0.63 
7 - 66 -0.52 


A OO l_ A 1 A 

0.2o ± 0.10 


A OOO 1 A Al a 

0.2o2 zb 0.01b 


91.5/97 




KaU J 05o2.y-oV01 


< r < a 


1 1 A 

1 — 10 


7 /IR+0.38 
7 - 46 -0.36 


A O A _1_ A H7 

O.o4 ± 0.0/ 


A 07yl _l_ A AAO 

0.2/4 ± O.OOo 


OOft 1 /OQA 

286.1/239 


A 07/) 7 

0.2/4/ 




0.5 < r < 4' 


0.4 — 10 


7 aq+0.31 
I -Uy_ .44 


A OA 1 A AO 

0.29 ± 0.08 


A OCO 1 A AAA 

0.262 ± 0.009 


OOft E /OAO 

226.5/202 




* 






c ao~t"0.35 

6.98I 028 


A OA 1 A AA 

0.30 ± 0.09 


A OCO 1 A A 1 1 

0.2d2 zb 0.011 


OO A O /OAO 

224.3/202 








1 — 10 


7 7C+0.72 


A OO 1 A 1 A 

0.28 ± 0.10 


A OTA 1 A A"l f* 

0.259 ± 0.016 


143. 7/137 




T~) v /~\ TAC £ O C E C T 

KACJ0o5o-5557 


< r < 8 


2 — 12 


i q cn+0.71 

13.59_ 05g 


A O A _l_ A AO 

0.24 ± 0.03 


A 007 1 A AAO 

0.287 ± 0.002 


COI O /CAO 

621.3/608 


0.2965 




0.5 < r < 4' 


2-12 


1456+2-21 

— I). 09 


0.21 ± 0.04 


0.290 ± 0.005 


378.1/400 




* 






14 61 +0 ' 57 


0.21 ±0.04 


0.291 ± 0.004 


377.1/400 




* 




0.4-10 


9 fjq+o.22 

t '- u ' :, -0.17 


0.23 ± 0.02 


0.288 ± 0.002 


1214.5/892 




RXCJ1131. 9-1955 


< r < 8' 


1 - 10 


ft 42 +o- 26 


0.26 ± 0.05 


0.299 ± 0.003 


322.5/283 


0.3075 




0.5 < r < 4' 


0.4- 10 


6 71 + 008 
"•'■ 1 -0.40 


0.23 ± 0.09 


0.285 ± 0.009 


263.1/255 




* 






6 RO+O' 38 
°- oz -0.21 


0.26 ±0.07 


0.287 ±0.012 


248.4/255 








1 - 10 


7.44«;^ 


0.26 ±0.07 


0.295 ± 0.024 


163.4/174 




RXCJ2337.6+0016 


< r < 8' 


1 - 10 


8.02«;f 7 


0.22 ±0.07 


0.328 ± 0.009 


211.1/255 


0.2779 




0.5 < r < 4' 


0.4 - 10 


6 77+ 29 
' ' -0.24 


0.24 ± 0.06 


0.315 ±0.011 


176.5/234 




* 






6.83t°; 2 2 2 8 


0.25 ± 0.06 


0.318 ±0.009 


172.7/234 








1 - 10 


7 EA+0.44 
'• du -0.37 


0.22 ±0.06 


0.317 ±0.011 


116.2/161 





'Introduce the upper limit of the residual Galactic emission described by an "apec" model in the double background subtraction 
method. 



able with the XMM-Newton observations of the REFLEX- 
DXL-like clusters. Similar analyses of the REFLEX-DXL 
clusters are in progress. RXCJ0307. 0—2840 is also very 
bright, but at a slightly lower redshift than the selection 
range for the REFLEX-DXL clusters. 

3.3.1. Gas distribution 

The regularity of the photon distribution shown in 
Fig. suggests that there are no large substructures in 
RXCJ0307. 0—2840. We thus assume a radially symmetric 
gas distribution. In order to get the actual gas distribu- 
tion we directly deproject the data from the spectroscopic 



analysis to get the gas mass profile (cf. Fabian et al 1981; 
Kriss et al. 1983). We divide the rings used for the temper- 
ature determinations into small subrings and fit the nor- 
malization for each subring separately, fixing the tempera- 
ture profile to the values obtained from the above spectral 
analysis, the metallicities to Z — 0.2Zq, the redshifts to 
z op t as given in Table ^ and the Galactic absorption to 
«h as given in Tabled We use a constant metallicity here 
since we did not detect a significant variation of metallic- 
ity with radius within the error limits of our analysis. In 
the soft band, the X-ray emission is almost independent 
of the temperature (Fabricant et al. 1980). The gas mass 
in each spherical shell is proportional to the square root 
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Table 5. Temperature profiles of REFLEX-DXL clusters obtained from the spectral analysis with residual background 
subtraction (double background subtraction) and without residual background subtraction. Col.(l): Cluster name. 
Col. (2): Energy band for spectral fit. Col. (3): Model used in XSPEC. Cols. (4-8): Temperature measurements. We 
do not obtain a consistent temperature measurement for RXCJ0528.9— 3927 in 4-8' region from the combined data 
because of the high background. 
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/cbJ l^e v j 










(keV) 




U < r < U.o 
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QflOO 
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mouen 
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i '- ot -1.07 
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b - 74 -0.29 


b - 5b -1.02 










moaelz 


o Q7+1-31 
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O.OO_0.55 
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R 4fi +2.14 
D- 4 D_ 1 4g 






1 


0-10 


model2 


10.4dIj 33 


10.69I 100 


n /IQ+0.65 
9 - 43 -0.59 


7 nn+0.47 
7 - 09 -0.43 


1 1 n/l+li-69 

11.94I 403 


i\AOJUU4o.4- 
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model2 


r oq+0.27 
O.OO_ 25 


a Qn+° 42 
b - 9U -0.37 


R 41 +0.32 
b - 4i -0.30 


r 90 + 0. 36 
u - z,:> -0.33 
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a mekal*wabs. 

6 mekal*wabs+powerlaw/b. 



of the integral emission, which can be calculated from the 
normalization of the spectroscopic analysis (e.g. Vikhlinin 
et al. 1999). 

In Fig. llUl we show the electron number density profiles 
for RXCJ0307.0-2840. Here, 1 arcmin = 0.240 Mpc at 
z = 0.2578. For the fits we use the standard (3 model 



n c (r) = n c0 



, r 
1 + I — 

Tc. 



(1) 



with the core radius r c and the shape parameter f3. The 
parameters from the best % 2 fits are listed in Table EJ The 
energy band used for the spectral fits does not affect the 
normalization, which corresponds to the electron number 
density profile (cf. Fig. I10[l . 



ture profile of RXCJ0307. 0—2840 drops towards the center 
due to cooling. We found that the parameterization 



k B T(r) 



1 



Ar 2 + Br + C 



(2) 



fits the measured temperature profiles quite well. Fig. ITT1 
presents the best x 2 fit for RXCJ0307.0-2840 with the 
parameters given in Table The temperature profile in 
the outermost regions can be fitted by a polytropic model, 
k B T(r) = fcBT (n c /n o0 ) 7 ~ 1 (e.g. Finoguenov et al. 2001b). 
In order for the system to be convectionally stable, the 
value of 7 should not exceed 5/3. Using the results of the 
spectral analysis in the 0.4-10 keV (1-10 keV) band we 
obtain 7 = 1.59 (7 = 1.46), which fulfills the stability 
criteria. 



3.3.2. Temperature distribution 

The precise estimate of the temperature structure greatly 
contributes to a reliable mass distribution. The tempera- 



3.3.3. Mass distribution 

We assume the intracluster gas to be in hydrostatic equi- 
librium with the underlying gravitational potential domi- 
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nated by the dark matter component. For the cosmological 
constant A = we have 



1 d{n c k B T) GM DM (r) 



pm p T 



dr 



(3) 



where n e and k B T are the electron number density and 
temperature distributions, respectively, and fi = 0.62 
is the mean molecular weight per hydrogen atom (e.g. 
Zakamska & Narayan 2003). 

Analytic models of the gas density and temperature 
profiles can be easily combined with Eq.© to obtain the 
mass profile: 



M(< r) 



k B T(r) 



3(3r 



(2Ar + B)k B T{r) 



•(4) 



We determine the mass distribution by using our tem- 
perature model and the /? gas density model. We use a 
Mont-Carlo simulation to calculate the error bars. 

Masses measured by the strong gravitational lensing 
are sometimes found to be larger compared to the mea- 
sured masses based on X-rays (Bohringer et al. 2000; Wu 
2000; Wu et al. 1998). High spatially resolved temper- 
ature profiles could help to resolve this discrepancy. In 
order to test the effects of temperature gradients we com- 
pare the mass estimates obtained under the assumptions 
of isothermality using the global temperature as measured 
in the 0.5 < r < 4' region fitted in the 1-10 keV band, 
and of non-isothermality. The mass and gas mass profiles 
are plotted in Fig.^] Under the assumption of hydrostatic 
equilibrium and isothermality, the virial radius and the to- 
tal gravitational cluster mass are 2.14 Mpc and 8.8 x 10 14 
M Q , respectively. 

3.3.4. Modeling a NFW mass distribution 

Navarro et al. (1997; NFW) described a universal density 
profile from numerical simulations in hierarchical cluster- 
ing scenarios, 



Pdm( 



^critPcrit 



(r/r s ){l + r/r s y 



(5) 



where S cr it and r a are the characteristic density and scale 
of the halo, respectively, and /? C rit, the critical density of 
the universe at the cosmic epoch z. <5 cr i t is related to the 
concentration parameter of a dark halo c = r V j r /r s by 



&-.rit — 



200 



3 ln(l + c) - c/(l + c) ' 



(6) 



We fit the observational temperature profile to obtain the 
parameters p s — <5 C ritPcrit and r s if we assume that the 
hot gas is in hydrostatic equilibrium with the dark matter. 
The former is well fitted by a standard (3 profile. The pa- 
rameters of the best fit of the NFW profile are presented 
in Table El The virial radius and virial mass estimates 



are smaller than the estimates under the assumption of 
isothermality. 

The NFW model describes the mass and gas mass frac- 
tion in the outer region well. Due to the cuspy NFW pro- 
file in the cluster center, the temperature fit based on the 
NFW model is higher than the observations. As a result, 
the gas mass fraction becomes lower in the center. But for 
this small central region we can not resolve the temper- 
ature structure well enough to perfectly recover the dark 
matter mass profile at the small radii. 

3.3.5. Modified hydrostatic equilibrium with A 

To be consistent with our background cosmological model 
with A ^ we should expect a second-order modification 
of the equation of hydrostatic equilibrium in the form 



1 d(n e k B T) GM DM {r) Ac 2 



/im p n c 



dr 



(7) 



The effect of a non-zero A is smaller than one percent 
and can thus be neglected compared to the relative error 
in our mass estimations. Sussman & Hernandez (2003) 
also point out a small effect of A on virialized structures, 
and that it could be significant only in the linear regime 
on large scales of r ~ 30 Mpc. 

3.3.6. Gas mass fraction distribution 

The distribution of the gas mass fraction is obtained ac- 
cording to the definition / gas (r) = M gas (r) / MdmO") • Since 
the derived gas mass is not completely unrelated to the de- 
rived total mass, we calculate the error bars of the gas frac- 
tion, A/ gas = ^(AAfgasAfoM) 2 + (M gas AM DM )7^DM- 

The profile of the gas mass fraction based on the NFW 
modeling is steeper in the central region as seen in Fig. 1131 
In both mass modelings the gas fractions increase with ra- 
dius, and range from / gas = 0.035 ± 0.012 to 0.138 ± 0.026 
in the outermost regions (r out = 1.441 Mpc) (see Tablc[SJ). 

These results are in good agreement with the measure- 
ments of Allen et al. (2002) based on Chandra observations 
of seven clusters yielding / gas ~ 0.105-0.138, and with the 
measurements of Sanderson et al. (2003) based on ASCA 
GIS & SIS and ROSAT/PSPC observations of 66 clusters 



yielding / g 



0.13 ±0.01 h 



-3/2 
70 • 



Our value is below the 



universal baryon fraction obtained with the recent WMAP 
measurement /b = fib/^m = 0.166, where ilb h 2 = 0.0224 
and fl m h 2 = 0.135 (Spergel et al. 2003). This reassures 
the estimate of the mass distribution. 

We also make use of the 8 < r < 10' region to obtain 
an upper limit of the gas mass in this shell applying the 
model "mekal+wabs" without residual background sub- 
traction. The upper limit of the total gas mass within 10 
arcmin (2.37 Mpc) is 2.49 x 1O 14 M . Based on the above 
mass modeling described by Eq. , the upper limit of the 
gas fraction within 10 arcmin is 0.49. This unreasonably 
high value confirms that the background dominates in this 
region. 
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Table 6. Parameters of each model from the best x 2 fits. r out = 1.441 Mpc is the outermost region where we can 
measure these parameters, rsoo and r v - n - are measured from the data. 



Model 


p SXcLlTlGt Gr 


0.4-10 keV 


1-10 keV 


Q 


1 C ^J-V-LJJl^ J 

n n ()Q~ 2 cm~ 3 ) 


0.061 ± 0.004 


0.062 ± 0.004 




2 43 + OS 


2 37 ± 04 




B 


60 + 02 


60 ± 02 


Eq.© 


A ( Mpc" 2 keV -1 ) 


0.283 ± 0.046 


0.331 ± 0.082 




B ( Mpc -1 keV -1 ) 


—0.240 ± 0.028 


—0 285 ± 052 




C ( keV -1 ) 


0.195 ± 0.004 


0.181 ± 0.007 


NFW 




267 + 028 


306 ± 034 




o. ho 15 MnMDC~ 3 1 


2.2 ± 0.3 


1.9 ± 0.5 




Q 


5 81 ± 59 


5 56 ± 68 




Msnn fl0 14 MrO 


4 10 ± 37 


5 13 ± 85 




renn fJVlDC) 


1.038 


1.116 




M„w fl0 14 MrO 


5 52 ± 49 


7.13 ± 1.18 




r vir (Mpc) 


1.550 


1.701 




/gas(^" < ?"vir) 


0.144 ±0.028 


0.124 ± 0.036 




M rout (10 14 M Q ) 


5.17 ±0.46 


6.30 ± 1.04 




/ga S (r < rout) 


0.138 ±0.026 


0.112 ± 0.032 


Eq.JU 


M 500 (1O 14 M ) 


6.86 ±0.01 


5.17 ±0.85 




r 50 o (Mpc) 


1.231 


1.118 




M ront (1O 14 M ) 


6.91 ±0.01 


16.46 ±4.69 




/gas(r < r out ) 


0.103 ±0.010 


0.035 ±0.012 



4. Summary and Conclusions 

We studied eight clusters of the REFLEX-DXL sample, 
selected from the REFLEX cluster survey at redshifts 
around z ~ 0.3, and one supplementary cluster at red- 
shift z = 0.2578. The data are from the MOS1, MOS2 
and pn detectors of XMM-Newton. The consistent results 
from the three detectors, obtained by excluding the ener- 
gies below 1 keV, give a good confidence in the applied 
method and provide tight constraints on the ICM param- 
eters like the temperature, metallicity, and redshift. 

Some of the clusters have been previously studied with 
ASCA, ROSAT and Chandra. The Chandra measurement 
of the temperature of RXCJ0658.5— 5556 includes mainly 
the central r < 3' region, where the temperature is high. 
We measure the global temperatures over a larger radial 
range 0.5 < r < 4' and r < 8', respectively. 

The parameter which best characterizes cluster mass, 
and which is most relevant for studies of the LSS and cos- 
mology, is the hot temperature of the bulk of the ICM. 
To avoid contamination by a possible central cooling core 
and by a possible soft excess or residual calibration un- 
certainties, we excluded the central r < 0.5' region and 
the softest part of the X-ray spectrum (< 1 keV) yielding 
reliable temperatures. These are the global temperatures 
that will be used in Paper I to derive the X-ray tempera- 
ture function for this sample. 

We obtained the spatially resolved X-ray temperature 
profiles for each cluster. For the determination of temper- 
ature profiles, the good statistics of the data allowed us 
to derive temperature values in five radial bins. In the in- 
ner regions an accuracy of better than 10-20 % can still 
be reached while the errors increase in the outer or out- 



ermost two bins. The temperature varies as a function 
of radius by factors of 1.5-2. The intracluster gas is cool 
in the center of RXCJ0232. 2-4420, RXCJ0307.0-2840, 
RXCJ0528.9-3927, and RXCJ1131.9-1955. No signifi- 
cant cooling gas is found with temperatures below 2 keV. 
In the outer region, the temperature drops at different 
levels. The differences of the temperature profiles in the 
cluster center may reveal that some clusters have relaxed 
cooling cores (but not all) and to some degree that we see 
the effect of non-gravitational processes. In this respect, 
it is remarkable that cooling cores are not only found in 
clusters with symmetric and regular X-ray images which 
might suggest a relaxed dynamical state, but also in the 
elongated, very disturbed cluster RXCJ1131.9— 1955. 

To study RXCJ0307.0-2840 in detail, wc find a model 
which fits the complex temperature profile of this cluster 
quite well. 

The mass distribution of this cluster, based on the pre- 
cise measurements of the distributions of the temperature 
and gas density, is similar to the mass distribution ob- 
tained under the assumption of isothermality within the 
region we can measure. We investigated the gas mass frac- 
tion of RXCJ0307. 0—2840 and found an increasing gas 
mass fraction as a function of radius, which is typical for 
most clusters. In the outermost region of the cluster, it 
is below the value of the universal baryon fraction. The 
uncertainty of the gas fraction is mainly caused by the 
temperature measurement. Therefore, a reliable determi- 
nation of the temperature profile is a key point to obtain 
the precise estimates of both the mass and the gas mass 
fraction. It plays an important role in the M-T scaling 
relation. 
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Fig. 2. Residual background (after subtraction of the back- 
ground obtained from the XMM-Newton observations of 
CDFS) of RXCJ0307.0-2840 for the 9.2-11.5' region from the 
pointing centers of the MOS1 (top), MOS2 (middle) and pn 
(bottom) scaled to 1 arcmin 2 area. The data are fitted by a 
power law model (solid lines). 
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Fig. 3. Temperature profiles for RXCJ0307.0-2840 with 
residual background subtraction (upper panel) and without 
residual background subtraction (lower panel) fitted in the 
0.4-10 keV band for MOS1 data (dotted lines), MOS2 data 
(dashed lines), pn data (dash-dotted lines) and combined data 
(solid lines). Additional solid lines connect the temperature 
measurements for the combined data. 
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Fig. 4. Comparison of the temperature esti- 
mates for RXCJ0014.3-3022, RXCJ0043.4-2037, 
RXCJ0232.2-4420, RXCJ0307.0-2840, RXCJ0528.9-3927, 
RXCJ0532.9-3701, RXCJ0658.5-5556, RXCJ1131.9-1955 
and RXCJ2337.6+0016 from MOS1, MOS2 and pn fitted in 
the 0.4-10 keV band with residual background subtraction 
(solid lines) and without residual background subtraction 
(dotted lines). 
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Fig. 5. XMM-Newton spectra of RXCJ0307. 0-2840 (top 
panel, fitted in the 1-10 keV band) and RXCJ0658.5-5556 
(bottom panel, fitted in the 2-12 keV band) extracted from 
the 0.5 < r < 4' region (pn in grey and MOS in black) with 
mekal model, considering the Galactic absorption and model- 
ing the residual background by a power law. The ratios of the 
observational data to the models are in the lower parts of the 
panels (offset zero for pn, +1 for MOS1, +2 for MOS2). 
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Fig. 7. The temperature measurements versus low-E cut-off 
from two different methods, the double background subtraction 
method (black) and the method applied in Arnaud et al. 2002 
(grey). 
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Fig. 8. Temperature profiles for RXCJ0307.0-2840 from spec- 
tral fits with (solid lines) and without (dashed lines) a geomet- 
rical deprojection. 




0.01 



E 



* 0.001 



0.0001 



n I I — I — I I I I I 1 1 — I I I I I I 



0.01 




_l I I I I'll 



Fig. 11. Measured temperatures of RXCJ0307.0-2840 fitted 
in the 0.4-10 keV band (dotted lines) and 1-10 keV band (solid 
lines). The corresponding curves show the best fits using the 
NFW profile (grey) and Eq.© (black). 
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Fig. 10. Measured electron number density profiles of 
RXCJ0307.0-2840 fitted in the 0.4-10 keV band (dotted lines) 
and 1-10 keV band (solid lines), respectively. The correspond- 
ing curves present the best fits using the standard j3 model 
with the confidence intervals (dashed curves). 




Fig. 12. Mass profiles (top) for RXCJ0307.0-2840 based on 
the temperature measurements fitted in the 0.4-10 keV band 
(dotted curves) and 1-10 keV band (solid curves) using the 
NFW model (grey) and Eq.@ (black). An additional dashed 
curve presents the mass profile under the assumption of isother- 
mality. The solid curve (bottom) with the confidence intervals 
(dashed curves) presents the gas mass distribution. 
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Fig. 13. Gas mass fraction distributions based on the tem- 
perature measurements fitted in the 0.4-10 keV band (dot- 
ted curves) and 1-10 keV band (solid curves) using the NFW 
model (grey) and Eq.@ (black). 



